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1. Introduction 

Large-A^ reduced models have been proposed as the non-perturbative definition of superstring 
theory. Especially, the IIB matrix model [||] is one of the most successful proposals. The IIB matrix 
model is formally obtained by the dimensional reduction of ten-dimensional ^ = 1 super-Yang- 
Mills theory to zero dimensions. In the IIB matrix model, spacetime is dynamically generated from 
the degrees of freedom of the bosonic matrices, despite the fact that it does not exist a priori in the 
model. Superstring theory is well-defined only in fen-dimensional spacefime, and if is an imporfanf 
quesfion how our four-dimensional spacefime dynamically emerges. Monfe Carlo sfudies of fhe 
IIB mafrix model have a possibilify fo shed lighf on fhis quesfion from a firsf principle calculation. 

The Euclidean version of fhe IIB mafrix model is obfained affer a Wick rofafion of fhe femporal 
direction. If has a manifesf SO(IO) rofafional symmetry which, if spontaneously broken, yields a 
spacetime compactified to lower dimensions. However, its numerical simulation has been hindered 
by the “complex action problem”, because the Pfaffian obtained after integrating out the fermions 
is complex in general. 

Apart from the matrix models of superstring theory, there are many interesting systems that 
are plagued by the “complex action problem”. Eattice gauge theories with a non-zero chemical 
potential are the ones that have attracted most of the attention in this context. In this work, we 
apply the “factorization method”, which was originally proposed in ref. 0 and generalized in 
ref. [^, to the Monte Carlo studies of the Euclidean version of the IIB matrix model. The IIB 
matrix model has also been studied analytically by the Gaussian Expansion Method (GEM) ii- 
Preliminary results of our Monte Carlo simulation are consistent with the GEM results and provide 
evidence that the factorization method is a successful approach to studying interesting systems that 
suffer from the complex action problem. 


2. Factorization method 


Generally, it is difficult to numerically simulate the complex action system 

Z = J 


(2.1) 


Since is not real positive, we cannot view it as a sampling probability in the Monte Carlo 

simulation. One way to calculate the vacuum expectation value (VEV) of an observable G is to 
use the reweighting {G) = Here, (• • •) and (• • • )o are the VEV’s for the original partition 

function Z and the phase-quenched partition function Zq = / dAe^^°, respectively. This is not an 
easy task since the phase E may fluctuate wildly. In order to compute {G) with given accuracy one 
needs configurations, where V is the system size. This is called the “sign problem” or 

the “complex action problem”. 

Yet another problem is that the important configurations are different for different partition 
functions. This is called the “overlap problem”. We are plagued with this overlap problem in 
trying to obtain the VEV (G) through the simulation of the phase-quenched partition Zq. 

The factorization method was proposed in order to reduce the overlap problem and achieve an 
importance sampling for the original partition function Z [Q We select the set of the observables 


Z = {G/ilk = 1,2,--- ,n}, 


( 2 . 2 ) 
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which are strongly correlated with the phase F. In the following, we define the normalized observ¬ 
ables ffic = ^kl{^k)o- We employ the factorization property of the density of states p(vi, • • • ,v„): 

,Xn) = {Y\5{Xk-ffk)) = ,Xn)w{Xi,--- ,X„). (2.3) 

k=l ^ 

The constant C = {e‘^)o is irrelevant in the following. p^^\xi,- ■ ■ ,v„) = (H^^i ^i-^k — ^k))o is the 
density of states in the phase-quenched model. w(vi, • • • v„) = {e‘^)x is the VEV in the constrained 
system 

Zx= [ dAe-^° f\5{xk-^k). (2.4) 

4 k=l 

When the system size V goes to infinity, the VEV’s are given by {ffk) = Xk, where (Jci, • • • ,jc„) is 
the position of the peak of p (vi, • • • , v„). This can be obtained by solving the saddle-point equation 

r ■ (2-5) 

V^— uXf^ ^ V J uXj^ V J 

When we properly choose the maximal set of the observables Z, we achieve effective importance 
sampling for the original partition function Z [^]. 


3. Euclidean version of the IIB matrix model 


We study the IIB matrix model |jT]], which is defined by the following partition function: 

Z = (3 2) 

where the bosonic part S'b and the fermionic part Sf are respectively 


5b = -^tr[A^,Av]^ (3.2) 

5f = ^tr(r/r„(^r^)„^[A^,i/rp]). (3.3) 

The bosons A^ (p = 1,2, • • • , 10) and the Majorana-Weyl spinors (a = 1,2, • • • , 16) are x 
traceless hermitian matrices. In the following, without loss of generality we set g^N = 1. The 
indices are contracted by the Euclidean metric after the Wick rotation. F^ are the 16 x 16 Gamma 
matrices after the Weyl projection, and ^ is the charge conjugation matrix. This model has the 
SO(IO) rotational symmetry. In ref. ||^, it is shown that the partition function is positive definite 
without cutoffs. 

This model is formally obtained by the dimensional reduction of ten-dimensional ^ = 1 super 
Yang-Mills theory to zero dimensions. The IIB matrix model has the yK = 2 supersymmetry 

S^^U^ = ie(^r^}y/, = SpU^=0, sPy/ = £- (3-4) 

For the linear combination and Sg^^ = i(dg^^ — Sg^'^), we have 

[ai“\5f]A^ = -2/5^''£(^F^)^, [5f\5fV = 0, (a,b=l,2). (3.5) 
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This leads to the interpretation of the eigenvalues of the bosonic matrices as the spacetime co¬ 
ordinates. Hence, the spontaneous symmetry breakdown (SSB) of the SO(IO) rotational symmetry 
is identified with the dynamical compactification of the extra dimensions. 

The order parameters of the SSB of the SO(IO) rotational symmetry are the eigenvalues 
{n = 1,2, • • • , 10) of the “moment of inertia tensor” 

Tpv = (3.6) 

which are ordered as Ai > A 2 > • • • > Aio before taking the expectation value. If (Ai), • • • , (Xd) grow 
and (Arf+i),--- ,(Aio) shrink in the large-A^ limit, this suggests the SSB of the SO(IO) rotational 
symmetry to SO{d) and hence the dynamical compactification of ten-dimensional spacetime to d 
dimensions. This scenario has been studied via GEM in ref. ijs]]. The results of the studies of the 
SO{d) symmetric vacua for 2 < r/ < 7 are summarized as follows: 

1. The extent of the shrunken directions r = limA^^oo (n = d + , 10) is ~ 0.155, 

which does not depend on d (universal compactification scale). 

2. The ten-dimensional volume of the Euclidean spacetime does not depend on d except d = 2 

(constant volume property). Eor the extent of the extended directions R = lim;v^~ VKi {n = 
1,2, - • • ,d), the volume is E = = /lo^ ^ 0.383. 

3. The free energy takes the minimum value at J = 3, which suggests the dynamical emergence 
of t/rree-dimensional spacetime. 

In ref. [Q], the six-dimensional version of the Euclidean IIB matrix model was studied via GEM, 
and the six-dimensional version also turns out to have these three properties. The same model was 
studied numerically in ref. ||7|], and the results are consistent with the GEM results. 

Next, we review the mechanism of the dynam¬ 
ical compactification of spacetime in the Euclidean 
IIB matrix model 0. Integrating out the fermions, 
we have 

J dy/e-^^ = Pf.^, (3.7) 

where = -ifabc{^ii)a^Ap is a 16(77^- 

1) X 16(77^ — 1) anti-symmetric matrix. The in¬ 
dices a,b,c run over 1,2,-•• ,77^ — 1, and fabc are 
the structure constants of SU(77). A^ are the co¬ 
efficients in the expansion A^ = with 

respect to the SU(77) generators . Under the 
transformation Aio —)• —A 10 , Pf./# becomes com¬ 
plex conjugate. We define the phase of the Pfaffian E as Pf.^ = |Pf.^|e'^. is real for the 

nine-dimensional configuration Aio = 0. When the configuration is rZ-dimensional {3 <d < 9), we 
find = 0 for m = 1,2, •• • ,9 — d, because the configuration is at most nine-dimensional 

up to the (9 — r/)-th order of the perturbations. Thus, the phase of Pf.y# becomes more stationary 
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Figure 1: The VEV (A„)o with respect to 
the phase-quenched partition fucntion Zq up to 
77 = 32. 
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for the lower dimensions. The numerical results in ref. 0 also suggest that there is no SSB of the 
rotational symmetry in the phase-quenched model. We calculate {X„)o numerically, where (• • • )o is 
the VEV with respect to the phase-quenched partition function 

Z() = I (3.8) 

We use the Rational Hybrid Monte Carlo (RHMC) algorithm, whose details are presented in Ap¬ 
pendix A of ref. 0. The result in fig. |I| shows that (A„)o converge to ~ 0.383 at large N for all 
n = 1,2, • • • , 10. This suggests that there is no SSB of the SO(IO) rotational symmetry, and that the 
result is consistent with the constant volume property. 


4. Results 


The model suffers from a strong complex action problem, and we apply the factorization 
method to this system. It turns out to be sufficient to constrain only one eigenvalue namely the 
choice of the set £ in eq. should be £ = {A„}. This is because the larger eigenvalues do not 
affect much the fluctuation of the phase. This choice of £ is similar to that of the six-dimensional 
version of the IIB matrix model 0. When we constrain A„, the eigenvalues • • • Aio take 

the small value, which corresponds to the SO{d) symmetric vacuum, with n = d + 1. This leads us 
to simulate the partition function of the constrained system 


= J dAe-^^\Pf^\5{x-l„), 


(4.1) 


which is simulated via the RHMC algorithm. The ratio = A„/(A„)o corresponds to the square 
of the ratio of the extents of the extended and shrunken directions (r/Z)^, in the SO{d) vacua with 
n = d+l. The saddle-point equation (2.5) is now simplified as 

1 


= - ^^logw„(v), where 


fh^\x) = ^log(5(x-A„))o, Wnix) = {e'^)n^ = {cosr)^^ 
dx 


(4.2) 

(4.3) 


in fhe large-A limif. (• • • )n^ is fhe VEV of fhe partition function Z„ We have {e^^)n,x = (cosr)„ 
because under fhe fransformafion Aio —)• —Aio the Pfaffian Pf./# becomes complex conjugate while 
the bosonic action ( 0^ and the eigenvalues of the tensor ( |3.6[ ) are invariant. The solution of the 
saddle-point equation ( p~^ ) jc„ gives the VEV (A„) = x„ in the SO(d) vacuum with n = d +\. 
Solving this saddle-point equation amounts to finding fhe minimum of fhe free energy 


■^SOid) (■^) 


-^logp„(v), where p„(v) = (5(x-A„)), 


(4.4) 


in fhe SO{d) vacuum wifh n = d + 1. The GEM resulf suggesfs fhaf fhe free energy fakes fhe 
minimum for fhe SO(3) vacuum. In order fo reduce fhe CPU cosfs, we focus on fhe n = 3,4,5 
cases, which correspond fo fhe SO(2), SO(3), SO(4) vacua, respectively. 

In fig. ^ (EEET) we plof logw„(v) for n = 4 up fo V = 16, where we observe a good scaling 
behavior af small x 


1 


logw„(x) 


~ —Z^n- 


(4.5) 
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The coefficients a„ and b„ are obtained for each N, by fitting the data. Then, we extrapolate the co¬ 
efficients a„,bn and obtain the large-A^ limit, which corresponds to = limAr^+oc ^logw„(x). 

This is represented by the solid line in fig. § (LEFT). 

The function (x) has a scaling behavior around 0.4 < x < 1 

(x) ~ g„ (x), where (x) = c„ (x - 1) + (x - 1 )^. (4.6) 


Subtracting this effect in order to reduce finite-A^ effects, we plot (x) — for n = 4 in fig. § 
(RIGHT). We find that the results scale reasonably well up to = 24 in the small-x region x < 0.4. 
This implies the hard-core potential structure at small x. In the six-dimensional version of the IIB 
matrix model, this effect is absent in the one-loop approximation [^, but is observed in the full 
model without one-loop approximation ||^. The intersection of and —^^n{x) 

represents the solution of the saddle-point equation (O). Fig. ^ (RIGHT) shows that the solution 
x„ is close to ^ ~ = 0.404- • • for n = 4. For « = 3,5, too, we have obtained similar results, 

and the solution x„ is close to 0.404. This is consistent with the “universal compactification scale” 
property. 




Figure 2: (LEFT) logW 4 (x) against x^ at small x. (RIGHT) ^f^^\x) — up to A = 24. Its intersection 
with — ^'I>„(x) gives the solution X 4 . 


Next, we compare the free energy ( p~4| ) for the 
SO{d) vacuum. The free energy at x = x„ is 

^SO{d)= J_ dxj^fl^\x)-j^\ogWn{xn), (4.7) 

with n = d + \. Due to the scaling behavior (p^), 
the first term of the r.h.s of eq. (El vanishes at 
large N. Thus we compare ^logw„(x„). From 
fig- |3[ we see that the free energy =^so( 2 ) is much 
higher than ^so(3) ^so(4) around x ~ 0.4. It 

is still difficult to determine whether the SO(3) or 
the SO(4) vacuum is energetically favored. More 
analysis will be reported elsewhere. 
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Figure 3: ^logw„(x) forn = 3,4, A= 12,16 
and n = 5, A = 12. 
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5. Conclusion 


In this work, we have performed Monte Carlo simulations of the Euclidean version of the 
IIB matrix model using the factorization method, in order to study the dynamical compactification 
of the extra dimensions. The results turn out to be consistent with the GEM predictions. We 
have seen that in the phase-quenched model there is no SSB of the SO(IO) rotational symmetry, 
and that the volume of spacetime is consistent with the GEM results. The function f!P\x) has 
a hard-core potential structure, and as a result of that, the computed shrunken dimensions are 
found to be consistent with the GEM results. Also, we have succeeded in finding that the SO(2) 
vacuum is energetically disfavored, compared to the SO(3) or SO(4) vacuum. The results of the 
Eorentzian version of the IIB matrix model, where (3-1-1 )-dimensional spacetime is found to expand 
dynamically [10], and the scenario discussed in this work, suggest that the physical interpretation 
of the Euclidean IIB matrix model needs to be further investigated. 
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